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Abstract 

In the spin ladder compound BiCi^POg there exists a decisive dynamics of spin excitations that 
we classify and characterize using inelastic light scattering. We observe low-energy singlets and a 
broad triplon continuum extending from 36 cm -1 to 700 cm -1 in (aa), (bb), and (cc) light scattering 
polarizations. Though isolated spin ladder physics can roughly account for the observed excitations 
at high energies, frustration and interladder interactions need to be considered to fully describe 
the spectral distribution and scattering selection rules at low and intermediate energies. More 
significantly, an interladder singlet bound mode at 24 cm -1 , lying below the continuum, shows 
its largest scattering intensity in interladder (ab) polarization. In contrast, two intraladder bound 
states at 62 cm -1 and 108 cm -1 with energies comparable to the continuum are observed with 
light polarization along the leg (bb) and the rung (cc). We attribute the rich spectrum of singlet 
bound modes to a melting of a dimer crystal. Our study provides evidence for a Z2 quantum phase 
transition from a dimer to a resonating valence bond state driven by singlet fluctuations. 
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Quantum spin ladders provide exceptional model systems, which enable us to explore 
rich one-dimensional physics, novel states of matter, and complex phase diagrams. lr3 The 
two-leg ladder is described by the leg (Jy) and rung (J±) couplings. For r = J±/J\\ 7^ 0, the 
ground state is given by a short-range resonating valence bond (RVB) state and elementary 
excitations are either dressed triplet excitations (triplons) or pairs of bound spinons with a 
finite gap.- 

Spin ladders have shown a variety of fascinating phenomena and turned out to have a 
relevance to diverse fields. In an applied field, they realize Tomonaga Luttinger liquids 
characterized by fractional S = 1/2 spinon excitations.- - — Weakly-coupled spin ladders 
demonstrate confinements of spinon excitations, which are a condensed-matter analogue of 
the confinement of quarks in particle physics.— The spin ladders show macroscopic quan- 
tum coherence of entangled spin pairs, holding a promise for quantum computation.—^ In 
hole-doped spin ladders, d-wave superconductivity has been reported and they provide a 
conceptual framework to understand unconventional superconductors.— 

Until now, isolated spin-ladder physics has been in the focus of research. To our best 
knowledge, the possibility of a quantum phase transition has not yet been experimentally 
addressed due to a lack of relevant materials. The spin ladder geometry may lead to two 
different routes to such transitions: Interladder couplings can induce a transition from a RVB 
state to a 2D Neel ordered state with classical spin- wave excitations, 15 while modifications 
of a leg to a frustrated J\ — J2 chain can tune a ground state as well. When J± = and 
J 2 > 0.24117Ji, two decoupled, dimerized chains are recovered. Upon turning on the rung 
coupling, RVB fluctuations are enhanced and eventually the dimer order melts into the RVB 
phase. 16 BiCu2PC>6 is an enticing material to address the crossover of a J\ — J2 chain to a 
spin ladder since it realizes a two-leg ladder with frustrating couplings along the legs.— 

The spin topology of BiC^POg is sketched in Fig. 1(a). The crystal structure of 
BiCu2PC>6 consists of complex ribbons with the basic unit being a folded pair of edge- 
sharing CuC>4 plaquettes. The neighboring plaquette pairs are connected by shared corners 
along the b axis while the next-nearest-neighbor plaquettes are linked by PO4 tetrahedra. 
The leg coupling Jy = J\ appears between the corner-sharing plaquettes while the rung 
coupling of J±_ = J4 involves the plaquettes of neighboring ribbons. The frustrating next- 
nearest-neighbor couplings J 2 and J' 2 along the legs run between the copper plaquettes, 
joined by another plaquette via a P0 4 tetrahedron. Finally, the sizable interladder coupling 

2 



J3 is formed between the edge-sharing plaquettes. The static magnetic susceptibility, the 
magnetic specific heat, the Knight shift, and the magnetization evidence a non-magnetic sin- 
glet ground state with a spin gap of A = 32 K.— _ — Geometrical frustration and anisotropic 
exchange and interladder interactions are invoked to explain the field-induced transitions 
between different magnetic phases, the reduction of the spin gap, and the impurity- induced 
magnetic ordering.— i22i2i In particular, the field-induced incommensurate magnetic ordering 
is beyond simple spin ladder physics. However, there exists no experimental knowledge about 
the possibility of a Z2 transition from a columnar dimer phase (with a translational sym- 
metry) to a RVB one (without a translational symmetry). A hallmark of such a transition 
is the presence of softening singlet modes. 35 Raman spectroscopy is a unique experimental 
method which can probe magnetic excitations in the singlet sector, information which is 
essential to understand the role of low-energy singlet fluctuations in driving a transition 
between two different singlet states.— 

In this paper we present magnetic Raman scattering of BiC^POg as a function of temper- 
ature and polarization. The most distinct feature is the presence of singlet bound states of 
two triplons in both interladder and intraladder polarizations. The intraladder singlet bound 
modes are the anticipated melting of a dimer phase through singlet fluctuations while the 
interladder singlet bound mode heralds the transition to a 3D magnetically ordered state. 
This demonstrates that BiCu2P06 is in the vicinity of two quantum critical points while 
being at the border line to dimer and 3D ordered phases. 

Single crystals of BiC^POg were grown by floating-zone method. The dimensions of 
the investigated single crystal is 4 x 2 x 1 mm 3 . The sample was kept in an evacuated 
optical closed-cycle cryostat with a temperature range from room temperature down to 3 K. 
Raman scattering measurements were carried out in quasi-backscattering geometry with the 
excitation line A = 532 nm of a Nd:yttrium aluminum garnet (YAG) solid-state laser. The 
incident power of 3 mW was focused onto the crystal surface with a spot diameter of about 
100 /im to avoid heating effects. To fully suppress Rayleigh and stray light scattering, the 
low-energy cutoff of Raman spectra was set to u = 10 cm -1 . Raman spectra were collected 
via a DILOR-XY triple spectrometer and a nitrogen-cooled charge-coupled device detector. 

Figure 2(a) shows Raman spectra at 3 and 280 K obtained in (aa), (bb), and (cc) po- 
larization, respectively. We identify 56 out of 60 Raman-active phonon modes as predicted 
from the factor group analysis {Prima space group); T Raman = 18 A g (aa, bb, cc) + 12B lg (ab) + 
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FIG. 1. (Color online) (a) Schematic of a frustrated spin ladder with the nearest neighbor leg 
coupling Ji, the rung coupling J4, and the frustrating next-nearest-neighbor leg couplings J2 and 
J' 2 ■ J3 denotes the interladder coupling, (b) and (c) Cartoons of two-triplets Raman scattering in a 
RVB picture, (b) a snapshot of the ground state; (c) locally excited singlet state of two neighboring 
triplet excitations depicted as a pair of the opposite arrows. 

18B2 9 (bc) + V2Bz g (ca). Upon cooling from 280 K to 30 K the phonon modes show a mod- 
erate hardening due to lattice anharmonicity. For temperatures below A = 32 K, a slight 
softening of the phonon frequency by 0.2-0.4 cm -1 is observed (not shown here). Since such 
spin-phonon coupling is a general feature of low- dimensional antiferromagnets,— we omit a 
further discussion of the phonon modes. 

Hereafter, our analysis focuses on magnetic excitations. For spin ladders, inelastic scat- 
tering of singlets by light occurs as simultaneous excitations of two neighboring singlets into 
a higher singlet state of two bound triplets since a Raman scattering process conserves a total 
spin of ^-component (S 1 * * = 0). This process is sketched in Figs. 1(b) and (c). Figure 1(b) 
illustrates a snapshot of a nonmagnetic spin singlet of two S=l/2 Cu 2+ spins. Low-energy 
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excitations are composite bosons of two-particle states, which are a linear combination of 
triplets {\t°), (t 1 ), l^ 1 }} = {| ti) + I it), I tt), I 11)}- The final state of Raman scattering 
thus corresponds to \t^)\t^ +1 ) + \t})\t^ +1 ) + )\t} +1 ). This two-triplet excitation is drawn 
as a pair of opposite arrows in Fig. 1(c). Here the arrows represent the excited triplets. If 
there is an attractive interaction between the excited triplets, a singlet bound state (tH" 1 
or t°t°) is formed in the similar way as mesons (qq) are created by the confinement of quark 
and antiquark. 

In the following, we introduce a polarization notation specific to a spin-ladder geometry. 
The a-,b-, and c-axis is perpendicular to the ladder plane, along the leg, and along the rung 
direction, respectively [see Fig. 1(a)]. Thus, (bb) [(cc)] polarization denotes the leg (rung) po- 
larization. Both (bb) and (cc) polarization are designated as intraladder polarization, (aa), 
(ab), and (ac) polarizations are collectively called interladder polarizations. The distinction 
between the intraladder and interladder polarization allows separating intrinsic spin-ladder 
magnetic behavior from 3D effects. 

In Fig. 2(a) we observe a broad, structured continuum at T = 3 K in intraladder (bb) and 
(cc) polarizations. The continuum evolves to a quasi-elastic scattering at high temperatures. 
The temperature and polarization dependence of the continuum leads us to assign them to 
the above- discussed two-triplon excitations. 

To allow a detailed study of the two-triplon energy, spectral weight, and shape, the 
phonons are subtracted. The polarization dependence of the resulting magnetic spectra 
at 3 K is presented in Fig. 2(b). The broad magnetic continuum with a multiple-peak 
structure extends from 36 cm -1 to 700 cm -1 . 28 For the leg (bb) and rung (cc) polarizations 
the broad continuum and two resonance-like peaks at 64 (S2) and 108 cm -1 (S3) (i.e. located 
above the two-particle threshold 2A) have complex lineshape with comparable intensity. In 
the interladder (aa) polarization the continuum is about four times weaker than in the 
intraladder polarization. The ratio of the integrated scattering intensity is I^ 9 ^ : /( run s) ; 
IW = 1 : 0.93 : 0.24. Since the magnetic Raman scattering intensity is given by / oc 
\{i\Si- Sj\j)\ 2 , it scales in each polarization with the strength of exchange interactions. Thus, 
we can estimate the ratio of exchange constants using the relation j[ legS> : j| ru ™ 9 ^ : J^ l > = 
RiegV I( le d) : R rung \/ I( run 9) : RnV where R represents Raman scattering interactions.— 
Assuming that R is proportional to the square of hopping matrix elements,— we obtain 
jf eg ^ : J^ ungS> : J£ l > = 1 : 0.89 : 0.35. This puts an upper limit to the strength of the 
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FIG. 2. (a) Comparison of Raman spectra of BiCu2P06 between 3 K and 280 K in (cc), (bb), 
and (aa) polarization, (b) Magnetic Raman scattering obtained in (cc), (bb), and (aa) polarization 
(after subtracting all phonon modes). The black arrow denotes a midband singularity (MBS) and 
Sj(i=l-3) stands for singlet bound states, (c) Evolution of magnetic Raman spectra with increasing 
temperature in a 3D plot. Inset: Temperature dependence of the two-triplon Raman scattering 
intensity, (d) Temperature dependence of the respective magnetic Raman spectra in a 2D plot. 

interladder coupling due to the additional leg couplings, J 2 and J' 2) suggesting that the 
interladder coupling is substantial. Further, we note that the sharp peak at 24 cm" 1 (Si) 
lying below the continuum is most pronounced in the interladder (aa) polarization. The 
three modes have ratios of -Bf/A = 1.01,2.72,4.59. The polarization selection rule and 
energy scale indicate that the two S2 and S3 modes pertain to an intraladder physics while 
the Si mode is related to an interladder one. 

Theoretical studies of two-leg ladders have shown that two-magnon bound states are 
always present in both the singlet and the triplet channel.—"— Based on the exact diagonal- 
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ization of a realistic spin model for BiCu2P06, Tsirlin et al, also suggest the presence of a 
number of singlet bound states at the zone boundary.— We therefore assign the three modes 
at 24, 64, and 108 cm -1 to anticipated singlet bound states. Below we will corroborate this 
interpretation by analyzing the temperature and polarization dependence of them. Here we 
mention that the resonance-like behavior of the bound modes points to their singularity-like 
nature and is attributed to a strong attractive interaction of neighboring triplets due to 
in-leg frustration. 

We recall that the singlet bound state has been observed in the two-leg compound 
Cai4_ x La3;Cu2404i with J±/J\\ ~ 1 via phonon-assisted two-magnon absorption.— The two- 
magnon singlet bound state has a two-peak feature assigned to van- Hove singularities in the 
density of states of the singlet dispersion at the Brillouin zone boundary and at k « tt/2. 
The spectral weight and shape of the two-magnon-plus-phonon absorption spectrum bears a 
close resemblance to our two-triplon Raman spectra [compare Fig. 3 of Ref.— and Fig. 2(a)]. 
In spite of the apparent similarity, however, the number and details of the singlet bound 
states differ, pointing to a significance of frustration and enlarged exchange degrees of free- 
dom in BiCu2PC>6. In contrast to the optical conductivity, Raman scattering measurements 
probe magnetic excitations without involving a phonon. Thus, the three peaks (Si, S2, and 
S 3 ) should be interpreted as separate bound states. 

We now turn to the two-triplon continuum at higher energies. The Raman response of the 
two-leg ladder was calculated by using the continuous unitary transformation method.— ~— 
For a finite rung coupling, the main feature is the existence of a midband singularity(MBS), 
which splits the continuum into two bands. This is attributed to a combined effect of 
a van-Hove singularity in the triplet-dispersion near k = and the orthogonal relation 
between excited triplets induced by triplet-triplet interactions. This spin ladder feature is 
not visible in Cai 4 _. r La a; Cu2404i due to sizable four-spin cyclic exchange interactions.— 1 ^ 
On the contrary, the studied compound exhibits a MBS at 455 cm -1 , which separates the 
secondary band at about 510 cm -1 from the main band centered around 300 cm -1 [see the 
black arrow in Fig. 2(b)]. This confirms that the magnetic excitations at high energies are 
largely described by isolated ladder-like physics. We remark that the simple spin ladder 
model also reproduces the magnetic behavior upon chemical substitution.— 1 ^ However, a 
closer look unveils a fine structure with multiple maxima at 254, 366, and 418 cm -1 in the 
leg {bb) polarization and 110, 229, 310, 377, and 439 cm -1 in the rung (cc) polarization 
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[see the arrows in Fig. 2(b)]. The strong polarization dependence of the spectral shape 
and the presence of a number of peaks within the continuum are beyond the spin ladder 
physics. Kohama et al. reached essentially the same conclusion based on the field-induced 
anisotropic, complex phase diagram. 24 We ascribe the fine features of the continuum to high- 
energy bound modes and/or van-Hove singularities related to the frustrated spin chain. Here 
we summarize the multi-facetted magnetic behavior of BiCu2P06 according to an energy 
hierarchy. At high energies the magnetic behavior approaches that of an isolated spin ladder 
as demonstrated by the MBS. At intermediate energies the intraladder bound modes show 
up as frustrated ID effects while at low energies 3D interladder interactions lead to the 
interladder bound mode. 



1.5 
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FIG. 3. Representative magnetic continuum at 160 K and in (bb) polarization, which is fitted 
to a sum of one Lorentzian and six Gaussian profiles. The gray Lorentzian spectrum denotes 
the quasielastic (QE) scattering response. Inset: The intensity of the quasielastic scattering is 
normalized by T 2 to be compared with the magnetic specific heat (C m ). The solid line indicates 
the calculated zero-field specific heat. 

In Figs. 2(c) and (d) the detailed temperature dependence of the magnetic continuum in 
the leg polarization is plotted. With increasing temperature the low-energy singlet bound 
modes dissolve into a continuum of excitations and then into quasi-elastic scattering upon 
further heating. This is easily understood as a thermal destruction of bound states. We note 
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that essentially the same temperature dependence of the singlet bound states is observed in 
the frustrated, two-dimensional Shastry-Sutherland system SrCu2(BOs)2.— In this system 
the multiple bound states are due to the localization of triplet excitations on a strongly 
frustrated lattice. 

With increasing temperature the two-triplon continuum is systematically suppressed 
while the spectral weight shifts toward lower energies. The integrated intensity of the mag- 
netic continuum, corrected by the Bose-Einstein factor, is shown as a function of temperature 
in the inset of Fig. 2(c). Upon cooling, the intensity increases monotonically and shows no 
saturation towards lower temperatures. This is due to a thermal population of triplets as 
found in other spin-gapped systems.— Despite the apparent resemblance, there is some dis- 
crepancy. In case of the spin gap systems, the magnetic continuum is gradually suppressed 
without appreciable changes in its spectral shape. Therefore, the softening continuum ob- 
served in BiCu2POg might be due to the damping of high-energy bound modes which lie 
within the background continuum and appear as multiple maxima. 

We will now turn to the low-energy quasielastic (QE) scattering, which provides infor- 
mation about the critical spin dynamics. In low- dimensional antiferromagnets, the QE 
scattering is due to fluctuations of the magnetic energy density; 45 



where E(k,t) is the Fourier transform of E{r) = —(J2i>j Jij^i ' Sj$(. r ~ r «))- Applying the 
fluctuation-dissipation theorem in the hydrodynamic and high-temperature limit, Eq. (1) is 
simplified to a Lorentzian profile, 



where C m is the magnetic specific heat and D is the thermal diffusion constant. Eq. (2) 
allows deriving C m by exploiting the relation C m (T) oc Iq E (T)/T 2 . For a quantitative 
analysis we fit the magnetic scattering to a sum of Eq. (2) (solid line) and Gaussian profiles 
(dashed lines). The representative fit is displayed in Fig. 3. The gray area, described by 
the Lorenztian profile, corresponds to Iqe- In the inset of Fig. 3 we show the temperature 
dependence of the derived C m (T). It shows a Schottky-type maximum at about 25 K, which 
is characteristic for a spin-gap system. The whole temperature dependence is fitted to an 
expression of C m oc T~ 3//2 exp(— A/T) to extract a spin gap.— The fitting result is depicted 




(1) 




(2) 
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as the solid line with the the spin gap of A = 36 K. The obtained value is very close to the 
spin gap determined by thermodynamic measurements. 




T(K) T(K) 



FIG. 4. (a) Temperature dependence of Raman spectra in interladder (ab) polarization on a 
logarithmic intensity scale. Temperature dependence of (b) frequency, (c) linewidth, and (d) 
normalized intensity of the singlet bound state at 24 cm -1 in interladder (ab) polarization (full 
circles) and at 62 cm -1 in leg ibb) polarization (open squares). The solid lines are fits to the data, 
see the text for details. 

Figure 4(a) shows the temperature dependence of Raman spectra in interladder (ab) 
polarization on a logarithmic intensity scale. We identify a resonance- like singlet bound 
mode at E% = 24 cm -1 . The singlet bonding energy is estimated to e s = 2 A — i^f = 
23 cm -1 . The lowest singlet bound mode shows up in all polarizations except the leg (bb) 
polarization. The scattering intensity is most pronounced in (ab) and (aa) polarizations, that 
is, in interladder scattering geometry. Since in the exchange light scattering mechanism the 
Raman operator is identical to the spin Hamiltonian, the intensity of the bound states should 
be observed with polarizations parallel to the dominant exchange path, i.e., along the leg 
and the rung direction. Therefore, the polarization selection rule indicates a substantial 
interladder coupling and the off-diagonal scattering matrix element, indicative of strong 
anisotropic exchange interactions. Here we note that the higher-energy singlet modes, S2 
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and S3 are observed only in intraladder polarizations, confirming that they are part of an 
isolated spin ladder. The intraladder singlet modes are interpreted as a softening singlet 
mode, which results from the transition from a dimerized to a spin ladder phase. 25 The 
melting of the dimer crystal corroborates that BiC^POg has a RVB ground state while its 
spin dynamics is governed by singlet fluctuations. Therefore, BiC^POg lies in the vicinity 
to a dimer phase. This is further supported by the numerical calculations that show a 
crossover from a spin chain to a spin ladder— and provides a convincing explanation for 
the field-induced anisotropic commensurate-incommensurate transition.— In addition, the 
interladder bound mode can be taken as an indication to the proximity to a 3D long-range 
ordered state from the RVB one. Future pressure work is required to test this possibility. 

The temperature dependence of the frequency, linewidth, and normalized intensity of 
the interladder and intraladder bound modes is compared in Figs. 4(b), (c), and (d). With 
increasing temperature both bound modes undergo a large hardening by 10-12 cm -1 toward 
higher frequencies and broaden enormously. At the same time, their scattering intensity 
decreases rapidly. Since the singlet bound states are susceptible to scatterings on thermally 
excited triplons, the reduction of the bound energy, the damping, and the decrease of the 
scattering intensity with temperature are related to the density of triplet excitations. The 
parameters of the bound modes are fitted to a single function, 1 — Aexp(— As/fcgT) where 
A is a constant. The frequency shift and the damping are reasonably well described by 
the above equation with As = 30 — 35 K. However, sizable deviation with a much smaller 
As = 11 — 14Kis seen in the temperature dependence of the intensity. This might be due 
to a broad triplon band and/or proximity to a quantum critical point. 

In conclusion, we have presented a comprehensive Raman scattering study of the mag- 
netic excitation spectrum in BiC^POg, a frustrated spin-ladder compound with interladder 
couplings. We observe the presence of both interladder and intraladder singlet bound modes 
at low energies as well as of a broad continuum of two-triplet excitations. The intraladder 
singlet bound states are a consequence of the melting of a dimer crystal into a RVB state 
and thus prove a Z 2 transition between two different singlet states. The interladder singlet 
bound state indicates an instability toward another ordered phase. BiCu2P06 is a model 
system, which hosts two different kinds of quantum phase transitions in a single system. 
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